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Abstract 

At  present  the  use  of  hypercomplex  methods  is  pursued  by  a growing  number  of  mathem- 
aticians, physicists  and  engineers.  Quaternionic  and  Clifford  calculus  will  be  applied  on 
wide  classes  of  problems  in  very  different  fields  of  science.  We  explain  Maxwell  equations 
within  the  geometric  algebras  of  real  and  complex  quaternions.  The  connection  between 
Maxwell  equations  and  the  Dirac  equation  will  be  elaborated.  Using  the  Teodorescu 
transform  we  will  deduce  an  iteration  procedure  for  solving  weak  time-dependent,  Maxwell 
equations  in  isotropic  homogeneous  media.  Assuming  the  so-called  Drude-Born-Feodorov 
constitutive  laws  Maxwell  equations  in  chiral  media  were  deduced.  Full  time-dependent 
problems  will  be  reduced  to  the  consideration  of  Weyl  operators. 

1 Historical  oriented  introduction 

Classical  Maxwell  equations  were  discovered  in  the  second  half  of  the  nineteenth  century 
as  result  of  the  stormy  development  of  electromagnetic  research  in  that  time.  The  study 
of  these  equations  has  attracted  generations  of  physicists  and  mathematicians  but  some 
of  their  secrets  are  still  hidden. 

At  about  the  same  time,  also  new  algebraic  structures  were  invented.  W.R.  Hamilton 
discovered  in  1843  the  algebra  of  real  quaternions  as  a generalization  of  the  field  of 
complex  numbers.  Under  the  influence  of  H.  Grassman’s  extension  theory  and  Hamilton’s 
quaternions,  W.K.  Clifford  created  in  1978  a geometric  algebra,  which  is  nowadays  called 
Clifford  algebra.  Its  construction  starts  with  a basis  in  the  signed  Rn  = Rp,q  with  units 
e\,  ...,e„.  Assume  that  ef  = -1,  for  i = 1,  and  e J = 1,  for  j = 1,  as  well  as  the 
anticommutator  relation 

e i ctj  -f"  ejCi  — 0 

for  i j-  Together  with  eo  = 1 one  can  construct  a basis  in  the  2" -dimensional  standard 
Clifford  algebra  Clp,q.  Incidentally,  in  1954  C.  Chevalley  [5]  showed  that  each  Clifford 
number,  i.e.  each  element  of  ClP:q,  can  be  identified  with  an  antisymmetric  tensor. 

Let  us  go  back  to  the  electromagnetic  field  equations.  Already  J.  C.  Maxwell  [15] 
himself  and  W.  R.  Hamilton  [10]  used  these  new  algebraic  techniques  to  try  to  simplify 
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Maxwell’s  equations.  The  aim  was  to  obtain  an  equation  of  the  type 

Du  + au  — F 

with  suitable  operators  D and  a.  For  this  reason  Hamilton  introduced  his  “N’abla  oper- 
ator” as  well  as  the  notion  “vector” . The  tendency  of  algebraisation  of  physics  continued 
in  the  first  half  of  the  last  century.  A long  list  of  important  publications  were  devoted 
to  this  topic.  We  only  stress  here  some  of  the  milestones,  beginning  with  the  “Theory 
of  Relativity”  by  L.  Silberstein  (1914)  [18]  , and  H.  Weyl’s  book  “Raum-Zeit-Materie”  of 
1921.  Important  results  of  Einstein/Mayer,  Lanczos  and  Proca  followed.  In  1935  this  de- 
velopment highlighted  with  the  thesis  of  M.  Mercier  (Geneva)  [16].  After  the  reinvention 
of  the  concept  of  “spinors”,  firstly  appeared  in  1911  in  a paper  by  E.  Cartan,  D.  Hestenes 
[11,  12,  13],  F.  Bolinder  [3]  and  M.  Riesz  [17]  wrote  fundamental  algebra  papers  with 
applications  in  electromagnetic  theory,  using  the  framework  of  Clifford  numbers  and 
spinor  spaces. 

Meanwhile,  in  the  late  thirties  the  famous  Swiss  mathematician  R.  Fueter  and  his  co- 
workers and  followers  used  a function-theoretic  approach  for  the  same  problems.  These 
ideas  were  refreshed  and  fruitful  extented  by  R.  Delanghe  and  his  group  and  A.  Sudbury 
in  the  seventies  and  early  eighties  (cf.  [4,  20]).  Influenced  by  the  success  of  complex 
analysis  and  Vekua  theory  a generalized  operator  theory  with  corresponding  singular 
integral  operators  [19]  and  a corresponding  hypercomplex  theory  for  boundary  value 
problems  of  elliptic  partial  differential  equations  were  developed  [8], [9]. 

Making  use  of  a transformation  of  Maxwell’s  equations  into  a system  of  homogeneous 
coordinates  we  will  propose  an  alternative  solution  method. 

2 Maxwell  equations 

Let  G be  a bounded  domain  with  sufficient  smooth  boundary  F that  is  filled  out  with 
an  isotropic  homogeneous  material. 

Using  Gauss  units  Maxwell  equations  read  as  follows: 


c rot  H 

= 4tt  J + dtD 

(Biot-Savart-Ampere’s  law) 

c rot  E 

= -dtB 

(Faraday’s  law) 

div  D 

= 47 rp 

(Coulomb’s  law) 

div  B 

= 0 

(no  free  magnetic  charge) 

Furthermore,  the  continuity  condition  has  to  be  fulfilled: 

div  J = -dtp, 

where  E = E(t,x)  is  the  electric  field,  H = H(t,x ) the  magnetic  field,  J = J(t,x)  the 
electric  current  density,  D = D(t,x ) the  electric  flux  density,  B = B(t,x)  the  magnetic 
flux  density,  p = p(t,  x)  the  charge  density,  and  c is  the  speed  of  light  in  a vacuum. 

The  relations  between  flux  densities  and  the  electric  and  magnetic  fields  depend  on 
the  material.  It  is  well-known  that  for  instance  all  organic  materials  contain  carbon  and 
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realize  in  this  way  some  kind  of  optical  activity.  Therefore,  Lord  Kelvin  introduced  the 
notion  of  the  chirality  measure  of  a medium.  This  coefficient  expresses  the  optical  activity 
of  the  underlying  material.  The  correspondent  constitutive  laws  are  the  following: 

D = e E + e /?  rot  E (Drude-Born-Feodorov  laws), 

B = p H + p (3  rot  H, 

where  e = e{t,  x ) is  the  electric  permittivity,  p = p(t,  x)  is  the  magnetic  permeab- 
ility and  the  coefficient  /3  describes  the  chirality  measure  of  the  material.  In  isotropic 
cases  one  has  the  possibility  to  use  the  so-called  Tellegen  representation 

D = e E + a H, 

B = pH  + a*E. 

The  connection  between  the  electric  field  E and  current  density  J is  given  by 
, J — o E + a g 

where  a is  the  electric  conductivity  and  g a given  electric  source. 

Starting  with  /3  = 0 and  replacing  D and  B by  D = e E and  B = pH  we  get  in  the 


case  of 

e = e{x)  , p — p(x) 

-edtE  + c rot  H = 47 rJ,  (2.1) 

pdtH  + c rot  E = 0 , (2.2) 

e div  E = AttP  - (Ve  • E)  , (2.3) 

pdivH  = -(Vp-H).  (2.4) 

After  summing  (2.1)  and  (2.4)  as  well  as  (2.2)  and  (2.3)  we  obtain 

—edtE  + c rot  H + p div  H = -(Vp  ■ H)  + An  J,  (2.5) 

pdtH  + c rot  E + £ div  E = — (Ve  ■ H)  + 4n  p.  (2.6) 


In  the  case  of  e,  p being  constants  we  can  introduce  the  new  functions  E.  H which 
are  defined  on  a homogeneous  space  with  a first  coordinate  ,Tq  and  the  other  coordinates 
x ='(®r,®2>*3)-  We  obtain: 

E(t,x)  =: 

The  equations  (2.5)-(2.6)  transform  into 

diE  + rot  H + p c div  H = An  J, 
diH  + rot  E + £ c div  E = An  p. 
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3 Quaternionic  representations 

Let  ei,e2,e3  be  the  generating  units  of  the  algebra  of  real  quaternions  TH . which  fulfil 
the  conditions 

a%ej  4-  e^e*  — 2<5 $j  (h  j — 1)2,3). 

This  leads  to  the  following  multiplication  rule  for  two  quaternions  u = uo+u  , v = Vo+V- 
UV  = UoVo  — u ■ V + UqV  + VoU  + U X V (Vi  £ Ft), 

where  u — uxex  + u2e2  + u3e3  , v = vxex  + v2e2  + 77363.  Further,  let  u = u0  + u be  a 
quaternion.  Then  u = uq  — u is  called  to  be  its  conjugate  quaternion.  The  operator 
defined  by 

D = diei  + d2e2  + d3e3 

is  called  Dirac  operator.  It  acts  on  a quaternionic  valued  function  as  follows: 

Du  = — div  u + rot  u + grad  uq  . 

With  the  multiplication  operator  m$ 

mou  — duo  + u (6  e !R+), 

with  u = uq  + u,  u = u\ei  + u2e2  + 77363,  we  obtain 

m^ddiE  + DH ) = 47 r J , 
mec{diH  + DE ) = 47 t p , 

and  so 

d\E  + DH  = J , 

diH  + DE  = m~c&Trp. 

Finally,  we  get 

d(E  + H)  = di(E  + H)  + D(E  + H)  = + m~}p)  =:  Fx  , 

d(E  — H)  = dx(E  - H)  - D(E  - H)  = 4t m^p)  =:  F2  , 

where  d is  also  called  Weyl  operator  and  d is  the  conjugate  to  d.  By  the  way,  a function 
u is  called  quaternionic  regular  if  d u = 0 and  quaternionic  anti-regular  if  d u = 0. 

For  simplifying  we  set:  E + H =:  v and  E — H =:  w.  Then  it  follows 

dw  = Fx(v,w),  (3.1) 

d w = F2(v,w).  (3.2) 

Let  us  have  a closer  look  at  the  functions  FX,F2.  The  electric  current  density  J is  given 
by 


J = a E + a g, 
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where  E and  g are  vector  functions.  This  leads  to  the  following  simplification 


Fi  = 47T 

F2  = Are 


<r(E  + g)  + — 

EC 


a(E  + erg) 


EC  J 


= 2w  |rr(v  + w)  + o g + 
2n  \a(v  + w)+g  — 


£ 

eel 

£ 

EC 


Hence 

Thus 


F2  — —F]  . 

dw  — Fi(v,w), 
dw  = —Fi(v,w). 


(3.3) 

(3.4) 


4 Integral  representation 

Let  G be  a bounded  domain  in  R3  and  a a positive  constant  . We  consider  in  IR  4 the  cyl- 
inder Z — Gx  [—a,  a) . A right  inverse  to  the  Weyl  operator  is  the  following  Teodorescu 
transform: 

(Tzu)(x)  = — / e(x  - y)u(y)dy  , Z = Gx  [-a,  a] 

<73  J 

z 

with  e(x)  = x/\x\i)  <r3  = 27r3/2/r (3/2) . We  obtain  in  a straightforward  manner 


dTzu 


in 

in 


Z, 

ra\z, 


and 


Tz  d u + <f>z 


in 

in 


R4  \ Z, 

with  <f>z  € ker  d.  In  complete  analogy  a conjugate  Teodorescu  transform  is  introduced. 
We  just  have  to  replace  e(x)  by  its  conjugate.  Now  it  follows  from  (7)-(8)  that 

v — TzF\(v,w)  + <t>z  {d(j)Z  = 0), 

1 w = TzF2(v,  w)  + 4>*z  {d(j)*z  = 0). 

Furthermore  we  have  to  introduce  Cauc.hy-Bizadse-type  operators,  which  are  defined 
by  the  boundary  data.  These  operators  read  as  follows: 


{Fozu){x) 


— f e(x  - y)n(y)u(y)d(dZ)y  (x  £ dZ) 

&3  J 

dZ 


and 


(Fdzu)(x)  :=  — [ e(x  - y)n{y)u{y)d(dZ)v,  (x  <£  dZ) 

. <73  J 

dZ 

where  n(y)  = (no  +n)(y)  denotes  the  unit  vector  of  the  outer  normal  on  dZ  at  the  point 

y- 
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It  can  be  proved  that 

4>*z  = Fazw  and  <l>z  = Fdzv  in  Z . 

It  should  be  noted  that  we  do  not  need  the  whole  trace  of  the  functions  w and  v 

on  the  boundary.  We  just  have  to  consider  these  parts  of  trzv  ( trzw ) which  are  lying 

in  the  corresponding  Hardy  space  of  functions,  which  permit  a quaternionic  regular 
(quaternionic  anti-regular)  extension  into  Z,  accordingly.  We  get  the  integral  equations 

v = 4-7TCT  Tz{v  + w)  + 47r  Tz(ag  + — ) + h,  (4.1) 

w = Ana  Tz(v  + w)  + 47r  Tz(ag  - — ) + h*,  (4.2) 

where 


h = Fgztrgzv  and  h*  = Fgztrdzvj. 

If  h,  h*  are  known  then  under  smallness  conditions  the  iteration  procedure: 

vn  = 4-7R TTz(vn-i  + wn-i)  + 4nTz(ag  + — ) + h, 

sc 

wn  = 4ircrTz(vn-i  + wn-i)  + 47rT|(f7c;  - £)  + h*, 

with  (vo  = Wo  = 0)  will  converge  in  suitable  Banach  spaces. 

Remark  4.1  In  [1]  is  proved  the  following  estimation: 

WTzWuL-co.C)  < 

5 Weak  time  dependent  Maxwell-equations 

Assume  now  e = s{x),p  = n(x),  k = k(x)  ( g = 0)  and 

E(t,  x)  = E0(t)E(x)  and  H(t,x)  = H0{t)Hi(x), 

where  the  scalar  functions  E0  and  H0  are  known.  Maxwell  equations  then  transform  to 


c Eq  rot  Ei 

= -dt(pH0)Hi, 

(5.1) 

c Ho  rot  Hi 

= ( dt{eE0 ) + 4ttkE0)E1  , 

(5.2) 

Eo(Ve  • Ei)  £ div  Ei 

= Anp, 

(5.3) 

(Vp  ■ Hi)  + p div  Hi 

= 0. 

(5.4) 

rot  Ei 
rot  Hi 
— div  Ei 


p dtH0 
c E0 
s dtE0 
c H0 


Hi  (XpHi 

47 TK  Ep\ 
c Ho) 


Ei  =:  (Ip Ei  , 


47rp  Ve 

FT  H 

eE0  £ 


= P 


-a- Ei  , 


It  follows 
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— div  Hi  = = . 

P - 

Here  a = ao  + a , 0 = 0o+  0,  or  :=  — ~ , (3  Using  the  fact  that  in  iff 

Du  = — div  v 4 rot  u, 

we  get 

D Ei  = oiqH\  -\-  p'  — ot-  E\ , 

D Hi  = PoEi-0-Hi. 

The  right  inverse  of  D is  the  corresponding  Teodorescu  transform  To  over  G C IR3.  A 
short  calculation  leads  to 

E\=  Tgo/oHi  — Tea  ■ E\  + Tap'  + 4>i, 

Hi  = Tg/3oEi  — Tg0  ■ Hi  + fa  , 
where  fa  € ker  D (i  = 1,2).  The  iteration  method 

E[n)  = -TGa  • E[n~1)  + TGa0H[n-1]  + Tap'  + fa, 

H[n)  = Tg0o  ■ E[n)  - Ta0  • H^~l)  + fa, 

with  H[0^  — f?[°^  = 0 converges  in  suitable  Banach  spaces  (L2)W2  ,C)  under  smallness 
conditions. 


In  the  time-harmonic  case  i.e.  Ho  = E0  = 1 and  e,  p,  are  constants  and  n = k(x)  we 
have 

DEi  = p'  and  DHi  — (30  E\ . 

Setting  0o  = 6~1  we  obtain 

D 6D  Hi  = ^ = p', 
i.e. 

A Hi  = -f. 

If  boundary  values  of  Hi  ( trrHi ) are  known  i.e.  trrH  1 = g the  complete  solution  is 
given  by 

Hi  = Frg  + TGVsDh  + TGQgdTGf  . (5.5) 

Here  Vs  and  Qg  are  orthoprojections  on  subspaces  in  the  quaternionic  Hilbert  space 
L2(G),  namely 

L2(G)  = 6 ker  D n L2(G)  0 D W2(G). 

s 
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The  scalar  product  is  defined  by 

( u , v)j  :=  J udvdG  € M. 

G 

The  operator  Vs  can  be  seen  as  a generalized  Bergman  projection. 

In  the  representation  formula  from  above  is  F?  the  Cauchy-Bizadse  operator  on  T 
and  h a smooth  continuation  of  g into  G.  Note  that  Vs  and  Qs  can  be  explicitly  defined 
(cf.  [9])!  Then 

E1  = ~ VsDh  + QsSTGf. 

47 XK 

Let  us  prove  that  the  boundary  condition  is  fulfilled!  Indeed, 

QsTgJ  = Df  with  / eW2  i-e-  trrf  = 0. 

TaDf  = / — Frf  = 0 (Borel-Pompeiu’s  formula). 

On  the  other  hand,  Plemelj-Sokhotzkij’s  formulae  yield: 

trrH\  = Prg  + trrVsDh  = Prg  + trrTDh  — trrTQsDh 
= Prg  + 9 ~ Pv9  + 0 = g. 

Pr  is  the  so-called  Plemelj-projection  onto  that  Hardy  space  of  IP-regular  extendible 
functions  into  G. 
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